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The post-recombination streaming of baryons through dark matter keeps baryons out of low mass 
(< IO^Mq) halos coherently on scales of a few comoving Mpc. It has been argued that this will 
have a large impact on the 21-cm signal before and after reionization, as it raises the minimum mass 
required to form ionizing sources. Using a semi-numerical code, we show that the impact of the 
baryon streaming effect on the 21-cm signal during reionization (redshifts z ^ 7-20) depends strongly 
on the cooling scenario assumed for star formation, and the corresponding virial temperature or mass 
at which stars form. For the canonical case of atomic hydrogen cooling at 10^ K, the minimum mass 
for star formation is well above the mass of halos that are affected by the baryon streaming and 
there are no major changes to existing predictions. For the case of molecular hydrogen cooling, we 
find that reionization is delayed by a /\z ~ 2 and that more relative power is found in large modes 
at a given ionization fraction. However, the delay in reionization is degenerate with astrophysical 
assumptions, such as the production rate of UV photons by stars. 



I. INTRODUCTION 

The relative velocity of baryons through cold dark mat- 
ter (v5c) after recombination is both supersonic and co- 
herent on scales of a few comoving Mpc [P-'T. This has 
potentially interesting implications for star formation [4- 

In particular, some baryons are able to stream out 
and escape the potential wells of mini-halos, reducing the 
available supply of gas to form stars [8^. Minihalos be- 
low 10^ Mq will have their baryon fractions strongly sup- 
pressed, and since the effect is velocity dependent (and 
therefore spatially dependent), there will be rare regions 
of higher streaming velocity where strong suppression can 
occur in halos as large as IO^Mq [8 . Overall, the es- 
cape of baryons from dark matter potential wells can be 
quantified simply by assuming a characteristic "filtering" 
mass, above which baryons can be captured in dark mat- 
ter halos [3 . 

Since the effect suppresses baryon capture by halos, 
one way to study it is through the 21-cm signal from hy- 
drogen during the Epoch of Reionization (EoR) . If stars 
reionized the universe, their formation history will shape 
the history of the ionization fraction and thus the 21-cm 
signal. As many observing teams are presently making 
rapid progress towards observing or constraining the 21- 
cm signal at the end of reionization [9^ ilOj , it is natural 
to ask what the imprint of the effect will be on this 
measurement. 

We find in this work that the result depends strongly 
on the assumed minimum mass for star formation. If the 
minimum mass needed to cool and produce stars is higher 
than the suppression scale, there should be (almost) no 
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effect on reionization. However, if the minimum mass 
needed to cool is below or of order the suppression scale, 
the vi)c effect will raise the minimum mass to the suppres- 
sion scale in a position-dependent way, and will couple 
the vi)c power spectrum to the 21-cm power spectrum. 
The traditional filtering mass for star formation is set by 
the linear theory evolution history, and thus it is domi- 
nated at low redshift {z < 150) by adiabatic cooling [llj . 
The virial temperature of a halo depends on its mass [12] , 
/i, and redshift z: 
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where /i is the mean molecular weight. Since virial tem- 
perature increases monotonically with halo mass, it is 
possible to convert a characteristic cooling temperature 
into a redshift dependent minimum mass that can be 
compared to the filtering mass: [4| 1121114] . 
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An important unanswered question is what form of 
gas cooling is responsible for producing the stellar pho- 
tons which reionized the universe. Molecular hydrogen 
cooling at 300 K should occur in primordial gas, lead- 
ing to the formation of metal- free Population III stars, 
but the Lyman- Werner (LW) UV background from those 
stars will also rapidly photo dissociate molecular hydro- 
gen throughout the universe [T5[ |T6] . In the absense of 
H2, atomic hydrogen cooling at 10^ K would be most 
relevant, and would set the minimum halo mass for star 
formation at M^in = 1O^M0[(1 + z)/10)]-^/2 [17 . 

Radiative feedback makes it difficult to know whether 
atomic or molecular cooling was dominant for the first 
stars and galaxies. The dynamic range of scales involved 
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is difficult to simulate, and there are many aspects of 
star formation in the early universe that remain unclear 
[18]. The relative importance of the cooling types may 
be spatially dependent [17 , leading in principle to a com- 
plex interplay of spatially varying minimum masses from 
both the effect and the spatially dependent LW back- 
ground. Reference |TT, for example, finds that the first 
halos which form through molecular cooling will destroy 
their own reserviors of H2. This makes atomic cooling 
more relevant to the minimum mass during the reioniza- 
tion epoch, so that the spatially- varying effect of molecu- 
lar hydrogen can be ignored. Some recent studies suggest 
that molecular hydrogen will dissociate quickly and that 
"self-shielding" in some regimes is even less important 
than previously thought [19) . 

These recent results support the notion that the atomic 
cooling scenario is more likely to be dominant during 
the reionization epoch than molecular cooling, which in 
turn will decrease the impact of the v^c effect on the 
21-cm signal. However, a considerable body of previous 
work has emphasized the importance of the v^c effect to 
the 21-cm signal using calculations based on molecular 
cooling [1, 2 . While it is valuable to consider the impact 
of the in its most extreme incarnation, we consider 
the likely impact of the v^c effect as a function of our 
uncertainty about the cooling scenario, so that the cases 
may be discriminated from each other. 





Atomic cooling at T^in=10''K 












Molecular cooling at T^in=300K 
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II. METHOD OF CALCULATION 

The streaming is coherent over a few comoving Mpc 
(cMpc) scales while star formation occurs at scales many 
orders of magnitude smaller. We study the imprint of the 
effect with a semi- numerical code 21cm FAST [13] with 
a 400 cMpc box with 4 cMpc resolution. Semi-numerical 
codes like 21cm FAST realize random cosmological ini- 
tial conditions, and use linear pertubation theory and 
analytic prescriptions for the relevant physics to evolve 
the volume over time. This volume can be made large 
enough to represent a fair sample of the universe. 

We modified 21cm FAST in order to incorporate the 
effect by self-consistently generating a coherent vi^c 
velocity in each 4 cMpc voxel. We then modify the min- 
imum mass locally in that voxel according to the pre- 
scription set by the filtering masses in Ref [3 . Our study 
is conducted in the regime where Tg ^ Tcmb^ which is 
valid for 2: < 20 [20^. This greatly simplifies our analysis 
because we do not have to consider the inhomogenous 
heating of the IGM by the first soures and the effect of 
the vijc on it. We assume the default cosmological pa- 
rameters for 21cm FAST with no halo finding. 



10* - 



1000 




10 15 20 



FIG. 2. The magnitude of the Vbc effect needed to modify the 
minimum mass of ionizing halos, giving a particular cooling 
temperature, as a function of redshift. The shaded regions 
correspond to la, 2cr, 3cr, and 4a fluctuations of the Vbc effect, 
from bottom to top. 



FIG. 1. The evolution of the minimum mass created by 
the Vbc effect, at one and four times the rms velocity of the 
Vbc (cr) (from Ref. 3 ), compared to the minimum mass for 
star formation associated with molecular and atomic hydro- 
gen cooling. 



We can then run realizations of the universe with the 
effect turned both on and off, for four different fiduci- 
ual cooling scenarios (and therefore minimum masses). 
The four scenarios involve minimum virial temperatures. 
Train, of: 200K, 300K, lO'^K, and lO^K. 

We obtain frequency and redshift dependent 21-cm dif- 
ferential brightness power spectrum, images, and mean 
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FIG. 3. Images of ionization fraction, 400Mpc on a side, with 
and without the Vhc effect. White represents neutral gas and 
black represents ionized gas. The first column has no Vhc ef- 
fect and is shown at a fixed Xi — 0.5, with each row using a 
different cooling temperature and redshift. The second col- 
umn uses the same cooling temperatures as the first column, 
but includes the Vhc effect and shows the equivalent panel at 
the same Xi — 0.5 (which for the top two rows is at a different 
redshift due to the delay in reionization) . The third column 
also includes the Vhc effect at the same cooling temperatures, 
but shows the image at the same redshift as the left column 
(which in the first two rows, is at a different ionization frac- 
tion). At higher cooling temperatures / minimum masses, 
there is essentially no difference in the three panels. 



ionization fractions as a function of redshift from 21- cm 
FAST. Our results are described below. 



modifies the image both at an equivalent redshift and at 
the same ionization fraction Xi = 0.5. As can be seen 
by comparing the left column images to the images in 
the right two columns, the v^c effect only imprints onto 
the ionization field if molecular cooling is important, in 
which case it increases the typical size of the neutral voids 
and delays reionization. Figure 4 shows the difference in 
power spectra due to the v^c effect, which clearly shows 
the transfer of power to longer wavelength modes. Figure 
5 shows the global delay in reionization due to the over- 
all higher effective minimum mass, in agreement with the 
results of Ref [2j . 
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FIG. 4. Ionization fraction history as a function of redshift, 
with and without the Vbc effect, with a minimum halo mass for 
star formation associated with cooling at 200K, 300K, lO^K, 
and lO^K (top left, top right, bottom left, and bottom right 
respectively) . The dashed line includes the Vbc effect while the 
solid line does not include the Vbc effect. Other parameters 
are 21cm FAST defaults (/* = 0.1, = 31.5). 



III. RESULTS 

We wish to find the redshifts at which the new filtering 
mass implied by the vi^c effect raises the minimum mass 
implied by the cooling scenario under consideration. In 
Figure 1, we show the redshift at which the minimum 
masses implied by molecular and atomic cooling become 
modified by a 1-sigma and 4-sigma fluctuation of the v^c- 
We also show in Figure 2 which regions of minimum-mass 
and cooling temperature parameter space are affected by 
the vi)c effect at different multiples of the rms velocity. 

To study the effect of the inhomogeneous v^c effect, we 
realize eight different cosmological histories using a mod- 
ified version of 21 -cm FAST, and compare the results as 
a function of minimum mass and by turning the v^c effect 
on and off. Figure 3 shows image slices of the ionization 
fields for the both cases, and compares how the v^c effect 



IV. CONCLUSION 

Given the likely scenario of atomic hydrogen cooling, 
the effect at z < 20 will probably not be detectable 
during the Epoch of Reionization. Molecular cooling may 
be dominant at z > 20, but in our fiducial model, sources 
which produce most of the UV photons causing reioniza- 
tion are formed at z < 20. It is very challenging to 
detect the corresponding 21-cm signal at z > 20, because 
the galactic synchrotron foreground has a brightness tem- 
perature T5 that rises steeply with observed frequencies 
Uobs, Tb (X 

The regime z > 20 also coincides with the time at 
which the first sources started heating the spin tempera- 
ture above the CMB temperature, making the 21-cm sig- 
nal both more interesting and harder to compute. How- 
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ever, the assumptions about cooling (and thus minimum 
mass) strongly modulate the vijc effect and future work 
should study how v^c predictions change in the high red- 
shift regime as a function of minimum mass prescriptions. 

If molecular cooling does turn out to be dominant for 
sources that produce most of the UV photons, the v^c 
effect will delay reionization. However, this delay will 
be degenerate with uncertainties in the star-forming ef- 
ficiency the escape fraction of ionizing photons from 
galaxies fesc^ and other astrophysical parameters such as 
the mass function of stars. We find that the shape of the 
power spectrum predicted by 21- cm FAST changes in a 
characteristic way from the no-Vbc case (at a fixed ion- 
ization fraction). Our ability to measure or constrain the 
effect in this case will depend on whether we can use 
simulations and semi-numerical models to reliably and 
accurately predict the 21-cm signal. 



FIG. 5. The 21-cm differential brightness temperature STb 
power spectrum, with and without the Vbc effect, at fixed 
x^=0.5, with a minimum halo mass for star formation associ- 
ated with cooling at 200K, 300K, lO^K, and lO^K (top left, 
top right, bottom left, and bottom right respectively). The 
dashed line includes the Vbc effect while the solid line does 
not include the Vbc effect. 
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